The Case for Hyperbolic Theories of Dissipation in Relativistic Fluids 
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I. INTRODUCTION 

As is well-known causal theories of dissipative fluids, 
both relativistic [Q and non- relativistic 0, || were set 
up to get rid of some undesirable effects connected to 
the conventional Eckart-type theories j3J. The advan- 
tages of the causal theories are the following ||]: (i) For 
stable fluid configurations the dissipative signals propa- 
gate causally, (ii) Unlike Eckart-type's theories, there is 
no generic short- wavelength secular instability in causal 
theories, (iii) Even for rotating fluids, the perturbations 
have a well-posed initial value problem. 

The central idea in causal-type theories (also called 
"hyperbolic" ) is to extend the space of variables of con- 
ventional theories by incorporating in it the dissipative 
quantities (i.e., heat fluxes, particle currents, shear and 
bulk stresses, etc.) of the phenomenon under considera- 
tion. Hence these quantities are treated on the same foot- 
ing as the "conserved" variables (energy density, particle 
numbers, etc.). The price to be paid is a more involved 
theory with a larger number of variables and parameters. 
In the approach of Israel [l.a-c] as well as in the one by 
Pavon et al [l.d], by combining the conservation of the 
stress-energy tensor and the conservation of the parti- 
cle current density with the Gibbs equation (which has 
extra terms due to the enlargement of the space of vari- 
ables) an expression is obtained for the entropy produc- 
tion. The transport equations for the dissipative fluxes 
follow from the latter by demanding that the aforesaid 
production be nonnegative. A distinguishing feature of 
these equations as compared with the corresponding ones 
in Eckart-type theories is that the former contain a "re- 
laxational term" -the convective time derivative of the 
corresponding flux multiplied by a parameter with the di- 
mension of time. Similar equations can be obtained in the 
extended thermodynamics approach of Muller and Rug- 
geri [2.b]. The transport equations can be understood as 



evolution equations for the dissipative variables as they 
(being "hyperbolic" in nature) describe how these fluxes 
evolve from an initial arbitrary state to a final steady 
one. On the contrary, the Eckart-type theories are fit for 
quasi-steady states only, and they are usually referred to 
as "parabolic" or "quasi-stationary theories" . 

The aforesaid time parameter r is usually interpreted 
as the relaxation time of the dissipative process. For in- 
stance, in heat conduction it is the time taken by the heat 
flux to attain its steady value since the establishment of 
a temperature gradient across the system. It reflects the 
fact that the effect (the heat flux or any other dissipative 
phenomenon) lags behind from its cause (the tempera- 
ture gradient or any other thermodynamic "force" act- 
ing on the system) . This is directly linked to the "second 
sound" -a well known phenomenon in laboratory physics 
whose justification can be found within hyperbolic theo- 
ries. We will turn to this topic below. Even in the steady 
regime the descriptions offered by causal and acausal the- 
ories do not necessarily coincide. The differences between 
them in such a situation arise from (i) the presence of r 
in terms that couple the vorticity to the heat flux and 
shear stresses. These may be large even in steady states 
(e.g. rotating stars). There are also other acceleration 
coupling terms to bulk and shear stresses and heat flux 
-see for instance [5. a]. The coefficients for these vanish 
in parabolic theories, and they could be large even in 
the steady state, (ii) From the convective part of the 
time derivative (which are not negligible in the presence 
of large spatial gradients), (iii) From modifications in 
the equations of state due to the presence of dissipative 
fluxes. However, it is precisely before the establishment 
of the steady regime that both types of theories (hyper- 
bolic and parabolic) can differ more importantly. There- 
fore if one wishes to study a dissipative process for times 
shorter than r, it is mandatory to resort to a hyperbolic 
theory which is a more accurate macroscopic approxima- 
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tion to the underlying kinetic description. Only for times 
longer than r it is permissible to go to a parabolic one 
(provided that the spatial gradients are not so large that 
the convective part of the time derivative does not be- 
come important, and that the fluxes and coupling terms 
remain safely small). All this has long been well-known. 

Recently, however, it has been argued that the physical 
content of causal (hyperbolic) and stationary (parabolic) 
theories is the same 0- Such a claim is based on the 
study of simple dissipative fluids under the assumption 
that its relaxation time is of the same order than the col- 
ision time between particles -something that, as we shall 
show is far from being always the case. Moreover in ]?J it 
is asserted that such a study is very general and should 
apply to essentially any theory of fluids (including those 
describing superfiuids, mixtures of different kinds of flu- 
ids, etc.); and that the meaningful differences between the 
causal theories and the non-causal Navier-Stokes theory 
can not be observed. The complicated dynamical struc- 
ture of the causal theories is necessary to insure that the 
fluid evolves in a causal and stable way. But this rich 
dynamical structure is unobservable, since the physical 
states of a fluid always evolve in a way that is also well 
described by the Navier-Stokes expressions for the stress- 
energy tensor, etc. These assertions (as well as some 
others of the like in || and 0) are misleading (to say 
the least) and likely to promt confusion. The purpose 
of |(| was to justify some class of relativistic, hyperbolic 
theories of dissipative fluids. This was achieved for sim- 
ple fluids under the assumptions that (i) the departure 
of these fluids from the Fourier-Navier-Stokes behavior 
is negligible small, and (ii) the relaxation times are ex- 
tremly short on macroscopic scales. Notice that in || it 
is not explictly denied that on examining a typical FNS 
(Fourier-Navier-Stokes) fluid in certain regimes, or with 
sufficient accuracy, there will be found non-FNS behav- 
ior, which could be satisfactorily described by some hy- 
perbolic theory. However, it is explicitly suggested that 
it may happen that the hyperbolic theories have no useful 
physical consequences at all. 

The aim of the present paper is to explicitly show 
that the physical content of hyperbolic (causal) transport 
theories is much broader than the corresponding one of 
parabolic (quasi-stationary) theories. We will do this by 
presenting a series of well-known physical results in an 
ample range of dissipative systems including simple flu- 
ids. A criticism in ||] and Q| of the hyperbolic theories 
is that they contain a number of arbitrary parameters. 
We would like to emphasize that this is not the case. 
In fact these parameters can be determined provided 
the state equations of the fluid under consideration are 
known. From these such parameters can be calculated 
either by kinetic theory or fluctuation theory. Israel and 
Stewart computed them for a quantum gas using Grad's 
14-moment approximation -see the set of equations (7.8) 
in [l.b]. For an example in radiative fluids, in which the 
coefficients connecting the heat flux with the gradient of 
the dissipative stresses are computed, see M. Further, 



a clear conceptual difference between the physical con- 
tent of both set of theories is the following one. With 
the help of Einstein's theory of thermodynamic fluctu- 
ations hyperbolic theories are capable of deriving the 
transport coefficients (heat conductivity k, bulk £ and 
shear n stresses, etc) from the fluctuations, because they 
are based on an extended entropy (which comprises the 
dissipative fluxes), at variance with parabolic theories. 
This has been done, for instance, for radiative fluids [||. 
This however cannot be done in the case of parabolic 
theories, which necessarily have to be supplemented with 
microscopic expressions to obtain these transport coeffi- 
cients. 

Section II recalls the concept of relaxation time. 
Section III brings about seven examples in the non- 
relativistic context showing that the physical content 
of hyperbolic theories is much wider than that of the 
parabolic ones. These include examples involving heat 
and/or electric conduction in solids. Note that there is 
nothing wrong in that because it is the electron and the 
phonon gases in the solid that are responsible for such 
transport processes within it JlO[ ] . Thereby it is meanin- 
ingful to apply, in this regard, fluid theories to solids. 
Section IV deals with four fully relativistic examples. In 
Section V our main conclusions and outlooks are sum- 
marized. 



II. THE RELAXATION TIME 

Before going any further it is expedient to say a few 
words about r. The relaxation time is a key quantity 
in hyperbolic systems with a distinct physical meaning 
as made explicit before. It is connected to the mean 
collisional time of the fluid particles t co i , sometimes erro- 
neously identified with it. In principle they are different 
since r is a macroscopic time, although in some instances 
it may correspond to a few t co i. It is more appropriate to 
interpret r as the time taken by the corresponding dissi- 
pative flux to relax to its steady value. We believe that 
at the root of the confusion raised by j(| and j?J is, among 
other things, the assumption that r and t co i are always of 
the same order. In fact one of the main arguments in Rcf 
H heavily rests on the assumption that the dimensionless 
quantity T = (tc s /L) 2 is always negligible. Likewise it is 
crucially needed for the bounds established by equations 
(3.22) and (3.23) in §] to hold. (Here c s stands for the 
isentropic speed of sound in the fluid under consideration 
and L the characteristic length of the system). That as- 
sumption would be right if r were always comparable to 
t co i and L always "large", but as we shall see there are 
important situations in which r 3> t co i, and L "small" 
although still large enough to justify a macroscopic de- 
scription. For tiny semiconductor pieces of about 10~ 4 
cm in size, used in common electronic devices submitted 
to high electric fields, the above dimensionless combina- 
tion (with r ~ I0~ 10 sec, c s ~ 10 7 cm/sec O]) can eas- 
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ily be of the order of unity. In ultrasound propagation as 
well as light-scattering experiments in gases and neutron- 
scattering in liquids the relevant length is no longer the 
system size, but the wavelenght A which is usually much 
smaller than L ||, p^ , G3. Because of this, hyper- 
bolic theories may bear some importance in the study of 
nanoparticles and quantum dots. 

Likewise in polymeric fluids relaxation times are re- 
lated to the internal configurational degres of freedom 
and so much longer than t co i (in fact they are in the 
range of the minutes), and c s ~ 10 5 cm/sec, hence it fol- 
lows T ~ 0(1). In the degenerate core of aged stars the 
thermal relaxation time can be as high as 1 second |L4j . 
Assuming the radius of the core of about 10~ 2 Rsun one 
has r ~ 0(1) again. Think for instance of some syrup 
fluid flowing under a shear stress, and imagine that the 
imposed shear is suddenly suppressed. This liquid will 
come to rest only after a much longer time (r) than 
the collision time between its constituent particles has 
elapsed. Many other examples could be added but we do 
not mean to be exhaustive. 

In transport theory when the mean free-path of the 
particles is much smaller than the macroscopic length 
scale of interest, the gas is in a state of small deviations 
from local thermal equilibrium (LTE) and the usual con- 
tinuum assumptions are appropriate. In the gas dynam- 
ical case the local thermal equilibrium equations are the 
FNS system which are parabolic in nature. When the 
mean free-path is of order of the macroscopic length scale 
the usual continuum description is no longer valid and 
the more fundamental kinetic approach must be consid- 
ered. Under rarefied situations (i.e., limiting to the one 
particle distribution function) the kinetic description is 
embodied in the Boltzmann Transport Equation (BTE). 
Solution techniques for this equation fall in two main 
groups: (i) direct, particle methods, and (ii) extended 
continuum approaches (e.g. based on moment methods). 
In general, direct particle methods (e.g. Monte Carlo) are 
capable of obtaining solutions under general conditions, 
whereas extended continuum approaches tend to be re- 
stricted to some definite ranges of the parameters. How- 
ever, for transitional fluxes ranging from very large mean 
free-paths (where only a direct particle method could 
be applicable) and very small mean free-path (where a 
usual continuum description is warranted) to utilize di- 
rect particle methods would lead to prohibitive compu- 
tational costs. Under these conditions extended contin- 
uum methods might be computationally much more ef- 
ficient. In the extended thermodynamics framework the 
moment equations obtained from the Boltzmann trans- 
port equation are closed by utilizing the entropy princi- 
ple (Muller et al [2.b]; Jou et al ||) or by assuming the 
maximum entropy distribution function (i.e., a special 
form for the distribution function which maximizes en- 
tropy subjected to the given moments ||). The resulting 
set of transport equations form a quasi-linear, symmetric 
and hyperbolic system of conservation laws (in a neigh- 
bourhood of LTE) which admits a supplementary con- 



servation law in the form of entropy balance. From the 
mathematical viewpoint this system has desirable prop- 
erties. The obtained mathematical model is able to de- 
scribe efficiently the regimes that lie in between the col- 
lisionless and the collisional regimes, i.e., the transitional 
one. In fact these models capture the usual continuum 
regime when the mean free-path is much smaller than the 
macroscopic length scale while in the transition regime 
provide values for the physical quantities that are con- 
sistent with a positive definite distribution function, and 
therefore are physically more satisfactory. These mod- 
els are potentially very useful for constructing hybrid 
continuum/kinetic schemes for simulation purposes (i.e., 
numerical schemes where the kinetic approach, which is 
computationally very expensive, is used only where it is 
required, that is, in the collisionless case, where contin- 
uum equations are used in the near LTE, and extended 
continuum models in the transition regions) p|. There- 
fore, hyperbolic theories broaden the domain of contin- 
uum models and bridge the gap between collisional and 
collisonless regimes. Then the scope of hyperbolic theo- 
ries is to describe the transitional regimes and it is mean- 
ingless to restrict them just to the collision dominated 
case. 



III. NON-RELATIVISTIC EXAMPLES 

1. Heat flux and shear stress in ideal gases 

Let us consider an ideal gas subjected to a temper- 
ature gradient and a shear stress. It has been es- 
tablished with the help of Grad's thirteen moment 
approximation [ fl6| that the the relaxation times 
for the heat flux and the shear pressure are given 
by T q — 3/(2nm7) and t s — 2r g /3, respectively. 
Here n stands for the particle number density, m 
the molecular mass, and 7 a positive constant re- 
lated to the interaction cross section -see §3.3 in 
Ref H . It is obvious that in a collision dominated 
regime n will be large and hence both relaxation 
times will not be much longer than t co \. In the 
opposite extreme (Knudsen regime) the continuum 
approximation breaks down. However, for interme- 
diate situations (where hyperbolic theories apply) 
both T q and t s can be substantially larger than t co l ■ 

2. Second sound in superfluids and solids 

The second sound was first discovered by Peshkov 
jl7j in Helium II in the mid-forties, and later con- 
firmed, also in solids, by a number of experimen- 
talists (see, e.g. |]|, f§] and @). It is related to 
the propagation of the heat signal in solids not at 
the phase speed predicted by Fourier's law, 

v p h = \/2xuj 

(which diverges as the frequency of the signal lo — > 
00), with 
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X = K/(pC v ) , 

where k, p and c v stand for the thermal heat con- 
ductivity, the mass density and the specific heat at 
constant volume of the solid, respectively, but at 
the speed 



Vph 



y/TL) + V 1 + ■ 



which remains finite at high frequencies. This re- 
sult straightforwardly follows from combining the 
Maxwell-Cattaneo equation, 
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typical of hyperbolic theories, with the equation for 
the balance of the energy in a rigid solid, 
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-see §6.1 of g| for details. At low frequencies (i.e., 
toj <C 1) the phase velocity of the hyperbolic theory 
reduces to the one predicted by the Fourier's law. 
However, outside that range both equations signif- 
icantly differ from one another; and in the high 
frequency limit (tlo 3> 1) the former reduces to 
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which is the propagation speed of thermal pulses - 
see §6.2 in ||. As mentioned above all these results 
have been confirmed experimentally -see |2(J and 
|H[ for a review. Notice that if one makes the un- 
physical assumption implied by equating t to zero 
(which corresponds to the Fourier's law) the latter 
expression diverges. 

Closely connected to the above discussion is the re- 
sponse of material sample to an instantaneous heat 
pulse (activated by means of a laser beam, or an 
electric discharge). An elegant treatment of the 
problem can be implemented by adding a delta-like 
term to the right hand side of the energy balance 
equation. This when combined with the Maxwell- 
Cattaneo equation leads to a second order differen- 
tial equation in terms of the associated Green func- 
tion, and one obtains for the temperature response 
to the heat pulse the expression 
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(where, for simplicity, we have restricted ourselves 
to one spatial dimension) while the temperature 



response corresponding to the parabolic theory 
(Fourier's law) is 



AT(x,t) 
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Here go is the energy supplied to the sample, £ = 
W( t / T ) 2 - ( x / Ut ) 2 , MO and J x (£) are related 
to the Bessel functions of the n kind by /„(£) = 
— (— i) n J„(i£), and H(U — t | x |) is the Heaviside 
unit step function. 

There are important differences between the predic- 
tions of the conventional parabolic theory and the 
causal theory, (i) In the former, the response to a 
heat pulse is instantaneously felt in all the sample. 
In the latter, by contrast, it vanishes for x > Ut. 
(ii) Two distinct contributions enter ATh yP erboiic- 
On the one hand there is the term with S (U t — x) 
which reflects the initial pulse damped by an expo- 
nential factor; on the other hand there is, in addi- 
tion, the Heaviside term which corresponds to the 
wake. The latter after a sufficiently long time re- 
duces to the conventional diffusion solution. Such 
a structure has been observed in the propagation 
of heat pulses in solids at low temperatures (see 
Figure 1 for a schematic representation of the ex- 
perimental results) and in the structure of X-ray 
peaks from X-ray bursters in star explosions -see 
point 1 of Section IV below, (iii) The hyperbolic de- 
scription predicts the presence of a precursor peak 
which has no counterpart in the parabolic theory. 
Although both situations are characterized by the 
same total energy, the distribution of the ther- 
mal energy differs from one another. While in the 
acausal description the energy is roughly evenly dis- 
tributed in the whole spatial region, in the causal 
description, by contrast, there is a sharp concen- 
tration of energy associated to the precursor peak. 

3. Ultrasond waves and light-scattering in ideal 
monoatomic gases 

In the low frequency regime (i.e., ™ <C 1) the 
sound propagation as well as light-scattering in 
ideal gases, with a heat flux and shear stress, are 
well described by the parabolic FNS theory of dissi- 
pative phenomena. However, as lo is increased the 
limitations of this theory begin to be felt (in fact 
as soon as tlo ~ 0.5 notorious discrepancies with 
respect to the FNS theory appear); and in the high 
frequency regime (™ 3> 1) the disagreements with 
experimental results become wider. Aside from 
predicting an unbounded sound speed as lo — > oo it 
fails to foretell the existence of an additional wave 
and yields erroneous values for the absorption co- 
efficient. Things fare differently for hyperbolic the- 
ories. These predict constant phase speed for both 
waves, 1.64 c s and 0.72 c s with c s = ^/5 fcsT/(3m), 
as lo — > oo p2|, pa]. All these results show good 
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agreement with the well-known experimental data 
of Meyer and Sessler and Greenspan (2j| about 
ultrasound waves. 

Table I contrasts the predictions of the FNS and 
hyperbolic theories |26| with the corresponding ex- 
perimental results in terms of the ratio between the 
mean free path I of the moleculas and the wave- 
length A of the soundwave. 

Light scattering in dilute ideal monoatomic gases 
has been studied in detail by Weiss and Miiller . 
These authors find good agreement of the hyper- 
bolic theory with the experimental data of Clark 
[^8| (see figures 7.1 through 7.3 for the dynamic 
form factor in ^7j) while the parabolic FNS theory 
fails hopelessly. 

4. Electric circuits 

Simpler and more familiar examples can be found 
in the province of electric circuits. To fix ideas con- 
sider a resistor and an inductance connected to a 
source £ (t). The intensity in the resulting circuit is 
governed by £{t) = RI{t) + Ldl(t)/dt. Note that 
if the source were suddenly turned off, the current 
will decay in a characteristic (relaxation) time given 
by r e = L/R. It is apparent that by suitably choos- 
ing R and L the aforesaid time can be made much 
longer than the collisional time of the electrons in 
the circuit. 

5. Thermohaline instability 

Imagine a layer of warm salt water above a layer 
of fresh cold water. This system will be dynami- 
cally stable against sinking so long the warmth of 
the salt water is high enough to reduce its specific 
weight to below that of the fresh water. However, 
as the upper layer cools off, its density augments 
and eventually small blobs of salty water will sink 
down. This kind of instability may also occur in 
stars [p9| . Let AT be the temperature difference 
between a generic blob of salt water and its sur- 
rounding and 
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the adjustment time where x denotes the mean ab- 
sorption coeffcient, p the mass density of the blob, 
a the radius assuming the blob spherical, as the 
Boltzmann radiation constant, and c p the specific 
heat at constant pressure. If the evolution of AT 
is analyzed using the Fourier law for the the con- 
ductive heat transport q = — kVT, one obtains an 
exponential decrease AT = AT(t = 0)exp(—t/Td), 
and the blob simply sinks in the water below it. 
However, if the Maxwell-Cattaneo equation is used 
in place of Fourier's one, it follows 

AT(t) = AT(t = 0)e- x \cos{ujx) 



2/3 \ sin(o;a;) 



Ar(0)r d 



with x = t/2r, lu = y^Ar/rd — 1 and 



AT{t)exp{t/r)dt 



and the blob undergoes an oscillatory motion -see 
[|30f for details. So for times shorter than the effec- 
tive time for relaxation into diffusion the behavior 
of the blobs differ drastically from that predicted 
by the parabolic theory. 

6. Polymeric fluids 

Simple models of dilute polymeric solutions ac- 
counting nicely for the observed dynamics of these 
fluids are, among others, those of Rouse and Zimm 
[pl| . In these the polymer is assimilated to a chain 
of beads connected by Hookean springs of constant 
K, and the viscous pressure tensor correspond- 
ing to the mode m is governed by the Maxwell- 
Cattaeno type equation 



P 



(m) 



i k 



Jm) Ur ik 

dt 



1 



dx k 



dvk 
dxi 



with ry*" 1 ) the viscosity coefficient and r'"^ = 
£/(2Ka( m }) the relaxation time of that mode. Here 
£ is the friction coefficient, a^ 71 ' = 4sin(7rr7j/2iV) 
the eigenvalues of the Rouse matrices and N the 
number of monomers. It is apparent that t*™) can 
be made as large as desired just by increasing this 
number. For large N one has r( m ) - N 2 . The 
above equation for Pju"^ can be easily derived in 
the framework of hyperbolic theories - see for in- 
stance Chapter 7 in Ref ||. 

7. Carrier transport in semiconductor devices 

Enhanced functional integration in modern elec- 
tron devices requires an accurate modelling of 
charge carrier transport in doped semiconductors 
in order to describe high field phenomena such as 
hot electrons, impact ionization, etc p^] . The stan- 
dard drift-diffusion equations |33| cannot describe 
adequately high field and submicron phenomena 
because they hold only in the limit of weak electric 
field and short mean free-paths The ultimate 
description (at least at the semiclassical level) is 
to find a direct solution to the semiclassical BTE 
PH by numerical techniques such as Monte Carlo 
simulation [fssf . Such an approach, however, re- 
quires very long computing times ad is therefore 
unsuitable for engineering design applications. An 
intermediate approach is to write a relatively small 
number of moment equations of the BTE, closing 
the system of equations by physical and mathemat- 
ical ansatz (the closure assumption) which must 
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then be checked by Monte Carlo simulation. The 
obtained macroscopic models, with the transport 
parameters (relaxation times, heat conduction co- 
efficient, etc.) obtained by Monte Carlo verified 
physical and mathematical ansatz, are called hy- 
drodynamical models and their range of validity is 
somehow in between the drift-diffusion regime and 
the kinetic one. 

Lately Anile et al [[36|, have introduced an ex- 
tended hydrodynamical model in which the closure 
of the moment equations is achieved by imposing 
the entropy principle of extended thermodynamics 
[2.b]. In the case of silicon semiconductors, there 
are several distinct time scales, the energy relax- 
ation time t w , the momentum relaxation time t p , 
the energy flux relaxation time r q , etc. Monte Carlo 
simulation shows that t w >> r p > r q , due to the 
predominance of elastic collisions, and this justifies 
the existence of a local thermal equilibrium state 
in which the electron temperature is much larger 
than the lattice temperature, and the electron gas 
average momentum and energy flux vanish. In the 
extended hydrodynamical model of Anile and Mus- 
cato pi , the first thirteen moment equations are 
considered and the closure, at the level of the flux 
of energy flux, is achieved by imposing compli- 
ance with the entropy principle up to second or- 
der around the state of local thermal equilibrium 
(supplemented by the maximum entropy ansatz for 
the distribution function in order to calculate one 
of the coefficients in the expansion, which other- 
wise would have remained undetermined |57| , ) . 
This closure has been checked by Monte Carlo sim- 
ulation in several cases. The resulting system of 
equations is hyperbolic in a broad domain of vari- 
ables. In particular the heat flux obeys an evolution 
equation of the Maxwell-Cattaneo type (and like- 
wise the viscous stresses, which however are unim- 
portant in this context and will be ignored). In 
the case of small gradients (of temperature, ve- 
locity, etc.) one can linearize the extended sys- 
tem around the state of local thermal equilibrium 
(the so-called Maxwellian iteration J39|), thereby 
obtaining a constitutive law for the heat flux . It 
turns out that the constitutive law for the heat flux 
comprises a convective term and a conduction one 
of the Fourier type. The resulting system of evo- 
lution equations for particle number, momentum 
and energy, closed with the above constitutive law 
for the heat flux, form a parabolic system which 
is consistent with the entropy principle of linear 
irreversible thermodynamics and the Onsager reci- 
procity relations. This system has been amply in- 
vestigated by various numerical techniques |40(| . 

This parabolic system should be contrasted with 
the extended hydrodynamical system of above, 
which is hyperbolic, and which comprises evolution 



equations for particle number, momentum, energy 
and an evolution equation for the heat flux. For 
both systems the relaxation times, which determine 
the transport coefficients and the production terms 
(arising from the collisions of electrons with the lat- 
tice phonons) , have been modeled in the same way 
and checked by Monte Carlo simulation in a ho- 
mogeneous slab of silicon. As a benchmark both 
systems have been solved numerically for the sub- 
micron diode consisting of a 0.1 /ito n + region 
followed by a 0.4 /im n region, and ending with a 
0.1 /im n + region. The doping density in the n + 
region is N = 5 x 10 7 cto~ 3 , while in the n region 
it is N — 2 x 10 15 cm~ 3 . The ambient device tem- 
perature is Tl = 300 Kelvin. We remark that the 
diode under consideration presents very steep den- 
sity gradients at the source and drain junctions and 
also very high electric fields (> 10 5 V/cm), which 
stretches up to its limits any continuum descrip- 
tion. Therefore this is an ideal benchmark for test- 
ing continuum models against kinetic simulations. 
The boundary conditions for the two models are 
consistent. The results for the velocity profile have 
been compared with Monte Carlo simulations with 
exactly the same conditions. It is apparent that 
the hyperbolic models, even in stationary situa- 
tion, perform much better than the corresponding 
parabolic one, vis a vis of the Monte Carlo simula- 
tions. Also for the velocity profile there is a spike 
in the parabolic model near the second junction, 
which is absent in the hyperbolic model (and also 
in the "true" Monte Carlo simulation -see Figure 2). 
The spurious unphysical result highlights the differ- 
ence between the parabolic and hyperbolic model 
in regions of high electric fields and large gradients. 

IV. RELATIVISTIC EXAMPLES 

Hyperbolic theories are conceptually more crucial in 
rclativistics scenarios as the "action at a distance" im- 
plied by parabolic theories finds no room there. 

1. Precursor peak in astrophysical X-ray bursters 
Many fluids of astrophysical interest fall in the cat- 
egory of " radiative plasmas" . By such it is under- 
stand an interacting mixture of a material gas (e.g. 
ionized hydrogen atoms plus electrons) and radia- 
tion quanta (e.g. photons). The material medium 
is considered to be locally in thermal equilibrium 
with itself as its internal mean free times are ex- 
tremely short, whereas the radiation quanta, be- 
ing out of equilibrium with the material medium 
(its mean free time is finite) gives to rise dissipa- 
tion. Obviously photons are scattered, absorbed 
and emitted by the gas ED. 

The analysis of elastic photon diffusion through a 
cloud of ionized plasma of modest optical depth via 



G 



parabolic theories predicts for the distribution of 
the photons over their time of escape from the cloud 
the conventional diffusion structure, but it is unable 
to account for the double peaked temporal luminos- 
ity profiles observed from X-ray bursters. Hoffman 
et al [^|| reported on a X-ray burster showing a 
precursor peak, that lasted about 4 seconds, neatly 
separated from the main event which lasted above 
1000 seconds. This pattern is reminiscent of those 
encountered in heat pulse experiments in solids al- 
luded before, it can be explained by the hyperbolic 
theory of photon transport [Q . The starting equa- 
tions are 



and 



N + J a =S 



1 



K T J a + J a + -h ab N, 



with N a (= Nu a + J a ) the photon four-current, 
u a the rest-frame four-velocity of the material 
medium, N — —N a u a the photon number density, 
ja _ faabjy^ ^ e con vective photon current rela- 
tive to the ionized plasma, h ab = g ab + u a u b the 
spatial projector, S the photon source (assumed 
Gaussian in time), and kt — ncrx where n is the 
electron number density and <tt the Thomson cross 
section. The solution of the corresponding disper- 
sion relation is given in terms of a Green function 
showing a diffusive behavior for times longer than 
t D r~j t * t co i/2 (the conventional diffusion time), 
and a well defined peak for times shorter than to- 
Here t* stands for the optical depth. As seen in 
figures 2, 3, 5, 6 and 7 of Ref E3] the solution is 
similar to that found in heat pulses in solids. 

2. Thermal relaxation in gravitational collapse 

Very often in the analysis of the gravitational col- 
lapse of fluid spheres the thermal relaxation time 
r is ignored 44 since it is usually very small as 
compared with the typical time scales of collapse 
for gravitating systems (for phonon-electron inter- 
action r ~ 10" 11 sec, and for phonon-phonon and 
free electron interaction r ~ 10" 13 sec at room tem- 
perature p5[). However, there are situations in 
which r cannot be neglected against the gravita- 
tional collapse. For instance, in cores of evolved 
stars the quantum cells of phase space are filled up 
and the electron mean free-path increases substan- 
tially and so does r. For a completely degenerate 
core star of radius one hundred of the solar radius 
and temperature of about ten million Kelvin one 
has r ~ 1 sec juj . In recent works Di Prisco et al 
J46| and Herrera and Martinez (47]] using the ther- 
mal conductivity by electrons in neutron star mat- 
ter and the Maxwell-Cattaneo transport equation 
for the heat flux, showed that the outcome of the 



collapse of the star is very sensitive to the value of r 
(they restricted themselves to the very conservative 
range 10~ 6 sec < r < 10 _3 sec), to the point that 
a fluid sphere may bounce if r is shorter than some 
critical value within the mentioned range, otherwise 
it undergoes a complete gravitational collapse -see 
figure 4 of Ref. |^7| . It is worthy of note that this 
analysis also applies (modulo some obvious adap- 
tation) to the Kelvin-Helmholtz phase of the birth 
of a neutron star whose duration is about tens of 
seconds and in which most of the binding energy 
(about 10 53 erg) is emitted Q. 

Further, phcnomenae prior to relaxation may influ- 
ence in a important way the subsequent evolution 
of the system (i.e., for times longer than r). In 
this connection it has been shown that for spheri- 
cally symmetric stars with a radial heat flow, the 
temperature gradient appearing as a result of per- 
turbations, and thereby the luminosity, are highly 
dependent on the product of the relaxation time by 
the period of oscillation of the star |^9| . 

3. Isotropic early Universe expasion 

It is natural to think that the fluid filling the 
primeval Universe was subject to steep gradients, 
whereby the very early cosmic world should be a 
right scenario for hyperbolic theories of dissipation 
to apply together with Einstein's field equations. 
As we will see by resorting to a couple of examples 
the predictions of the causal theories greatly differ 
from that of the acausal (parabolic) ones. Unfor- 
tunately there is no compelling observational evi- 
dence favoring one or another, as nowadays cosmic 
observation is still at its infancy. As a first example 
let us consider a homogeneous and isotropic uni- 
verse with geometry described by a flat Robertson- 
Walker metric with scale factor a(t). And assume 
that the source of the metric is a barotropic fluid 
(i.e., one with equation of state P — (7 — l)p , 7 = 
constant > 0), endowed with a nonvanishing bulk 
viscosity whose coefficient depends on the energy 
density via £ = ap m with a a positive semidefinite 
constant and to a constant parameter. By using 
in this setting the parabolic Eckart's theory of ir- 
reversible phenomena Barrow ]5Q] found the first 
order differential equation for the evolution of the 
Hubble parameter 



2H + 3(7 - 3 m aH 2m ~ 1 )H 2 



0, (H=-) 
a 



For to < 1/2 this equation predicts, among other 
things, that the Universe begins with a Friedmann 
singularity and approaches a inflationary expansion 
(i.e., H = constant) for t — > 00. However, if 
to > 1/2 one has what it is termed as a "defla- 
tionary" behavior, i.e., the Universe starts with an 
inflationary expansion and evolves toward a Fried- 
mann state. 
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Things are different when a hyperbolic theory is 
used instead. In such a case an expression for 
the relaxation time must be postulated. Pavon 
et al [|T) choose the well motivated relationship 
t == C,j p , later justified by Maartens |52]. The cor- 
responding differential equation for the evolution of 
the Hubble factor 



a(3 n H 2n H 



(1 + ial3 n 1 H n+m )H 
12na 



-H 



n-f m 



H* = , 



becomes second order in time, and so the num- 
ber of possible evolutions greatly augments. (Here 
n = m — 1, /3 is a short hand for 3/ (8ir), and units 
have been chosen so that c = G = 1). Last equa- 
tion reduces to the former one under the unphysical 
assumption of vanishing r, and in the very partic- 
ular case that H = —3-fHH. The set of solutions 
of the causal equation is far richer than that of the 
acausal equation, both for nearly stationary solu- 
tions and for the nonstationary ones -see figures 1-5 
in Ref fil] . Here we want to stress that the defla- 
tionary solutions found by Barrow |pp| for to > 1/2 
are seen to be unstable. Initial inflationary expan- 
sions can be found only by fine-tunining the ini- 
tial conditions H(0) and H(0), hence it is not a 
general feature of cosmic evolution. Another point 
of interest refers to the "generalized second law of 
thermodynamics" , which enters into play when the 
spacetime posses an event horizon. In such cases 
is not the entropy of the cosmic fluid that musn't 
decrease with time (Sf > 0), but the combined en- 
tropy of the event horizon Sh and that of the fluid, 
i.e., Sh + Sf > 0. The conceptual trouble with 
theories allowing unbounded speeds for dissipative 
signals is that no meaningful event horizon can be 
ascribed to the spacetime even though it experi- 
ences a de Sitter expansion. So one is forced also 
in this regard to resort to hyperbolic (causal) theo- 
ries if one wishes to deal with cosmic horizons when 
the source of the gravitational field is a dissipative 
fluid. The concrete expressions for the generalized 
second law, equation (3.53) of Ref and equation 
(53) of Ref ]|lj , largely differ from one another. 

4. Anisotropic early Universe expansion 

The above study has been extended to anisotropic 
spaces, Bianchi-type I |55| and Bianchi-type III 
[ p4f , which admit shear stresses. In both instances 
the dynamical equations turned out far too compli- 
cated to allow analytical integration, so qualitative 
analysis and numerical integration were carried out 
instead. For Bianchi-type I the initial anisotropy 
dies away quickly and, in general, neither the Fried- 
mann nor de Sitter expansions are stable. Never- 
theless, for a wide range of values of the param- 
eters occurring in these equations, there exists a 



stable submanifold of phase space where the latter 
presents asymptotic stability. This contrasts with 
the findings of Huang J5j|, who used Eckart's the- 
ory of dissipative processes [4. a] for the bulk and 
shear stresses, and concluded that the Bianchi-type 
I space asymptotically evolves either to Friedmann 
or de Sitter states. For Bianchi-type III the Fried- 
mann expansion is unstable and the de Sitter ex- 
pansion stable pM. 



V. CONCLUDING REMARKS 

By resorting to some specific examples in non- 
relativistic as well as in relativistic nonequilibrium ther- 
modynamic situations we have made explicit the fact that 
hyperbolic (causal) theories of transport have a range of 
applicability much broader than the parabolic (acausal) 
ones. The differences between both sets may be ignored 
in steady regimes only, i.e., when the relaxation times 
of the processes under consideration are much shorter 
than the time scale of interest. This has been known 
long since. It was originally hinted by Maxwell ]56|] , and 
first seen in the laboratory by Peshkov E^] and later 
confirmed many times over in variety of experimental 
situations |pl[ , H. In this way we have dispelled the 
claims recently made in j6| and j?J implying that hyper- 
bolic theories have no useful applications. Note that in 
any case the findings in |6) and Q are limited to simple 
fluids and situations such that the quantity T, defined 
above is much lower than unity. Thereby it does not 
cover many interesting physical fluid situations, such as 
those recalled in this article. Therefore claims such as hy- 
perbolic, relativistic theories for dissipative fluids, reduce, 
physically, to Navier-Stokes found in p[ cannot be justi- 
fied. In summary, the suggestion in [3| and [Q that the 
class of solutions to which their theorems apply include 
all solutions of the hyperbolic dissipative fluid theories 
having physical relevance, is at variance with experimen- 
tal evidence. 

Current hyperbolic theories are not free, however, of 
some drawbacks, which will be mentioned in the follow- 
ing. This does not imply that one should go back to the 
acausal parabolic theories, but instead try to improve 
on the existing hyperbolic models. The main drawback 
of the moment methods currently used in present day 
hyperbolic theories is the introduction of spurious char- 
acteristic speeds, i.e., characteristic waves that arise only 
from the truncation procedure (the number of moments 
retained in the theory), and have no counterpart in the 
underlying kinetic description. Since these characteris- 
tics constitute a fan which broadens as the number of mo- 
ments increases, there could be serious difficulties (due to 
the mathematical properties of the hyperbolic systems) 
of the appearance of unphysical discontinuities in the so- 
lutions of the hyperbolic theories at places where the 
fluid speed crosses one of the characteristic speeds Wf\. 
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Fortunately the problem is not so severe because lately 
Weiss [|58| has shown, by direct numerical integration (for 
the case of a rarefied gas of Maxwell molecules) that, by 
some algebraic cancellation, the discontinuities can oc- 
cur only at the highest characteristic speed, and the lat- 
ter increases with the number of retained moments. If 
this result would be confirmed (and mathematically un- 
derstood) in more general cases, the aforesaid difficulty 
would be alleviated. 

More generally, the root of the difficulty lies in using 
the usual polynomial moments. In fact an avenue for fur- 
ther research should be to investigate whether a different 
basis function would be more suitable. We remark that 
any basis function must be able to capture phenomena 
varying on a very short time scale. In fact, the require- 
ment of relativistic causality itself, forces the basis func- 
tion to be capable of describing high frequency situations. 
One formulation of relativistic causality is that the wave 
speeds must be less or equal to the light speed in vacuum 
c. This must be applied to the characteristic speeds, but 
also to the dispersion waves arising from a linearization 
of the equations. In this case the relativistic causality 
takes the form v g < c asw-xxj (where v g is the group 
velocity), and obviously this requires that the equations 
must be able to describe high frequency phenomena (at 
least in the asymptotic sense). This is the main physi- 
cal reason why quasi-stationary theories, based on slowly 
varying variables, are intrinsically unable to comply with 
relativistic causality and makes it mandatory to consider 
hyperbolic theories. 
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TABLE I 



l/X 


0.25 


0.50 


1.00 


2.00 


4.00 


7.00 


{Cs/ v ph)FNS 
{c s /v p h)Hyp. 
{Cs/VphjExp. 


0.40 
0.52 
0.51 


0.26 
0.43 
0.46 


0.19 
0.44 
0.50 


0.13 
0.47 
0.46 


0.10 
0.48 
0.46 


0.07 
0.49 
0.46 
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Table I 

Numerical values of c s /v p h as a function of l/X. 
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Figure 1 

Schematic results of heat pulse experiments. The dashed line represents the heat pulse at the source, whereas the 
solid line shows the structure of the heat signal at an arbitrarily selected point in in the sample. Peaks A, B and C 
correspond to ballistic phonons with longitudinal sound velocity, ballistic phonons with transverse sound speed, and 
second sound, respectively. The tail D corresponds to the diffusive heat propagation. This structure can be explained 
by hyperbolic theories. By contrast, parabolic theories are able to explain the tail D only. 



Figure 2 

Velocity profile for the electrons in a submicron diode under a very high electric field. The crosses correspond to the 
Monte Carlo simulation, the solid line to the prediction of the hyperbolic theory, and the dashed line to the parabolic 
one. The latter shows a spurious peak. 
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